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Abstract 

We correct the error in the sign of the pseudoscalar pole contri¬ 
bution to the muon g — 2, which dominates the 0(a 3 ) hadronic 
light-by-light scattering effect. The error originates from our over¬ 
sight of a feature of the algebraic manipulation program FORM 
which defines the e-tensor in such a way that it satisfies the rela¬ 
tion = 24, irrespective of space- 

time metric. To circumvent this problem, we replaced the prod¬ 
uct ^^^2(13^4 £viV2i , 3i'i by rjviv-i'O ij .2^2 V /o ^3 V i ■'' hi the FORM- 
formatted program, and obtained a positive value for the pseu¬ 
doscalar pole contribution, in agreement with the recent result ob¬ 
tained by Knecht et al. 
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In this brief article we report the result of our reexamination of the pseudoscalar 
pole contribution to the muon g — 2. In the previous studies |], |2|, the pseudoscalar 
pole contribution had been noted to be the dominant term of the 0(a 3 ) hadronic 
light-by-light scattering contribution to the muon g — 2, a p = (g p — 2)/2. In view 
of the expected accuracy of the new muon g — 2 measurement at the Brookhaven 
National Laboratory (BNL) (Aa M (exp) = 4 x 10~ 10 ), the present authors and 
Bijnens et al. have examined this contribution carefully [|, [3|, taking account of 
the experimental data obtained for the Py*y-vertex (P = 7T°, ?/,?/) at CLEO ['4J. 
The primary purpose of this paper is to examine the sign of this contribution in 
light of the two recent papers || ^J. We therefore concentrate on the 7r° pole 
contribution which gives the value 

^(tt 0 ) = -55.60 (3) x HT 11 , (1) 


in the naive vector meson dominance (nVMD) model |IJ. This model simply mod¬ 
ifies the Wess-Zumino term for the 7T°yy-vertex by attaching p-meson propagators 
which carry photon momenta /y, & 2 , 
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Both Bijnens et al. [|2|] and Bartos et al. [[7| obtained the negative value for the 
pseudoscalar pole contribution independently of our result. Taking this value 
into account as a part of the standard model prediction, it is found that the 
current measurement of g — 2 indicates 2.6 a discrepancy from the prediction of 
the standard model || Q. 

However, the recent papers II pointed out that the pseudoscalar pole contri¬ 
bution is positive, which is opposite to the sign of ([!]). If this is true, it will reduce 
the deviation of the current BNL measurement from the standard model consider¬ 
ably. In view of these papers !,|, and of the significance of the 0(a 3 ) hadronic 
light-by-light scattering contribution in interpreting the current measured value of 
a p , we decided to reinvestigate the 7r° pole contribution. 

In the following, we summarize the results of our investigation to show various 
points we have scrutinized closely. A full account is being prepared in a separate 
paper [ITU . 


The first phase of examination consists of the following four steps of the original 
calculation |1]|: 
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(1) We obtained the expression for the contribution to the /i/i 7 -vertex function 
F v {viiPf) (pi = p — 9/2 and pi? = p + q/2 are the initial and final momenta of the 
muon) by direct application of Feynman rules. 

(2) We evaluated the trace of the 7 matrices using the algebraic manipulation 
program FORM to extract a^ from T u (pi,Pf) by means of the magnetic moment 
projection [], 


S = i im n r n 21 Tr 9)F v (Pi,Pf)) , 

0 {p-Q—>0; p z -\-q z /4:-^m z } 

P u (p, q) = 16 ^ 4 q 2 [i> - | + rnj (( 7 - frY) p 1 2 ~ 3 q 2 p u ) 
xU + ^ + mj , 


( 3 ) 


where q is the incoming photon momentum. 

(3) The result of Eq. (^|) is plugged into a FORTRAN program written in the 
formalism developed for the numerical evaluation of Feynman diagrams [|ll[ . 


(4) The numerical evaluation of is carried out with the help of the Monte Carlo 
integration routine VEGAS [[ 12 ] . 


In addition we confirmed by hand calculation and by MATHEMATICA that 
the projection operator in Eq. (||) works correctly, by extracting the anomalous 
magnetic moment = Fm{ 0 ) from 

1 


r "(p!,p F ) = (Wb" + M9 2 )^W‘V). 


with 


^= 2^71 • 


( 4 ) 

( 5 ) 


After going through these steps 


(I) we rederived the value given in Eq. (P for the 7 r° pole contribution in the 
11 VMD model. 

The next phase was to derive the result obtained in Ref. ||. 

We first assumed that the expression given for in Eqs. (3.4) and (3.5) of 
Ref. || is correct and see if the numerical evaluation of them gives the value of 
our Eq. m but with opposite sign, as claimed in Ref. ||. We translated Eqs. (3.4) 

1 We correct the typo-error for the expression for P v (p, q) of eq. (3.24) in Ref. Q], which does 

not affect the result obtained there. 


3 







and (3.5) into an expression suitable for the numerical evaluation according to the 

By carrying out the five-dimensional integration over the 


formalism of Ref. 11 


Feynman parameters with use of VEGAS, we obtain the value in Eq. (|l|) but with 
an opposite sign. This test not only confirms the result of Ref. [5j but also gives 
an evidence that 


(II) both Ref. H and our work performed the loop integration part correctly. 
The disagreement in sign must therefore come from an earlier stage. 

We therefore switched our attention to the task of projecting out from 
F 1/ ( P i,Pf) by means of algebraic programs. Using the algebraic manipulation pro¬ 
gram FORM we pursued the steps described in Ref. || in detail to see whether 
we can reproduce their Eqs. (3.4) and (3.5). We found that 

(III) the trace operation by FORM yielded the sign opposite to that of 7j ) 2 (qi, q 2 ; p) 
in Eq. (3.5) of Ref. [|. 

We thus recognized the following two possibilities as most likely: 

(a) Ref. IQ] made a mistake in picking up the sign of the trace of the 7 matrices. 

(b) We made a mistake at the stage of picking up the sign of the trace of the 7 
matrices systematically. By systematically we mean that we failed to identify the 
sign irrespective of what projection operator was used to extract a /t ; we derived 
the value of Eq. (jl]) by both our own projection operator (|3]) and those given in 
Eqs. (2.9) - (2.11) of Ref. 0. All projectors gave the same result. 

I 11 examining the possibility (b), we noticed one crucial difference between 
Ref. 0, which leads to the positive value, and ours, which leads to the negative 
value; while Ref. 0 used the algebraic manipulation program REDUCE to perform 
the trace calculation of the 7 matrices, we used FORM instead. Recall that we 
have used FORM even for examining the results of Ref. p]. 

Thus we decided to check whether we handled FORM properly. The program 
FORM had been used successfully to calculate the QED corrections to the g — 2 
of the muon and the electron as well as other observables by one of the present 
authors. However, this does not guarantee that we deal correctly with the e-tensor, 
the central object of our study of the pseudoscalar contribution. 

A simple test of this question is to see if our naive use of the FORM declaration 
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(Fixlndex 1 : —1, 2 — 1,3: —1;) works successfully to verify the identity 

= - 24 , ( 6 ) 

which should hold in Minkowski space-time. Unfortunately, the result turned out 
to be +24. This means that this simple declaration does not work for the e-tensor 
in Minkowski space-time Q bi order to circumvent this difficulty we avoided the 
direct use of e-tensor and repeated the trace calculation using the right-hand side 
of the identity in Minkowski space-time: 

- ['Wi W 2 W 3 r W 4 ± ( tlie other 23 terms)] , (7) 

where the other 23 terms are obtained by shuffling the order of {z/ 1; z/ 2 , z/ 3 , i/ 4 } in 
all possible ways. Each term contributes in the bracket with the sign + (—) if the 
even (odd) permutation is performed to reach this order from {iq, i/ 2 , ^ 3 , ^ 4 }- We 
found that the method using ( 0 ) led to a result opposite in sign to the previous 
result ([[]). On the other hand, REDUCE passed the same test without difficulty. 
Hence, we conclude that it was the case (b) that actually happened. In other 
words, the 7 r° pole contribution in ([I]) must be changed 

a M (7r°) = 55.60 (3) x KT 11 , (8) 


in the 11 VMD model. 

The sign of the pseudoscalar pole contribution (||) has returned to that of 
Ref. |U| and agrees with that of Ref. ||. But, in the former case, the program 
containing the product of two e-tensors has also been used to extract a /t . At the 
same time, there was an error in the sign of the logarithmic term dominating 
the pseudoscalar pole contribution, which was corrected in Ref. |T| ^|. The double 
switching of the sign has led accidentally to the positive value for the pseudoscalar 


contribution in Ref. [13 


We conclude that the sign of all the results for the pseudoscalar pole contri¬ 
bution as well as the axial-vector meson pole contribution described in Ref. |], H 
must be reversed. The signs of the charged pseudoscalar loop contribution and the 
constituent quark loop contribution are not affected by the problem noted here. 
Collecting those results from Eq. (1.8) and changing the sign of Eqs. (1.9) and 

2 We thank J. Vermasseren for explaining how to use the e-tensor for the Minkowski metric. 

3 The delicate point will be described in Ref. [|Ioj . 
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(5.1) of Ref. we obtain the new value 

CV(LL) = 89.6 (15.4) x 10" 11 . (9) 

as the current value of the 0(a 3 ) hadronic light-by-light scattering contribution 
to the muon g — 2. This reduces the discrepancy between the measurement and 
the prediction of the standard model to 1.6 a deviations. 
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